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At  the  request  of  the  organizing  committee,  I 
Mould  Ilka  to  share  aooa  of  ay  observations  and 
reasabranoes  about  tha  developaent  of  computa¬ 
tional  aathaaatlea  In  Czechoslovakia  and  tha  USSR. 
My  obaarvationa  will  ba  vary  aubjaetlva  and  broad 
la  aoopa. 

A.  Tha  Development  la  Csaohoalovakla 
1.  Tha  davalopaant  until  1918 

A  vary  asaantlal  alias tona  in  tha  davalopaant 
of  aolanoa  In  Control  Europe  was  tha  foundation  of 
Uio  Charlao  Uni varsity  In  Prague  In  13*8.  To  ay 
knowladca,  tha  first  aathaaatloal  text  st  this 
Uni  varsity  mss  llksly  ‘HT-ffMS  Proaavcua  by 
KrlSfcan,  froa  Praehatloa  (InCxechoelovakla)  wrlt- 
tan  In  1A00.  Thla  taxt  oonoontrataa  on  arithmetic 
and  so  X  saa  It  as  tha  flrat  taxt  on  eoaputatlonal 
aathoaatlos  In  Control  Suropa. 

Many  outstandlnc  aathoaatlelans  Intar as tad  In 
computations  ware,  dlraotly  or  Indirectly,  for 
shortar  or  longer  parloda,  assoclatad  with  tha 
Charlas  Uni  varsity.  Let  aa  aantlon  tha  aatrono- 
aars  T.  Braha  (15*6-1601),  J.  Kaplar  (1571-1630) 
and  J.  Burgi  (1552-1632),  aaong  othars.  Tha  sil- 
var  alnlng  In  Bohaala  (tha  aajor  Mining  plaoa  In 
Europo  at  thla  tlaa)  and  tha  oonatruotlon  of  a 
aystaa  of  ponds  In  Southorn  Bohaala  required  alg- 
nlflcant  affort  and  high  aoouraey  In  goodoalo 
aaaauranants  and  ooaputatlons.  This,  togathor 
with  tha  naod  of  aatronoay,  contributed  to  tha 
davalopaant  of  eoaputatlonal  aathaaatlos.  Con- 
putatlonal  aathods  of  Braha  (how  to  Multiply  nun- 
bars  by  additions  with  halp  of  tablaa  of  aln  and 
eon)  togathor  with  tha  logarithms  (tables  of 
Maplar,  Kaplar,  Burgi),  and  tha  davalopaant  of  a 
aaohanloal  ooaputar  by  Sohlekart,  froa  Tubingen 
In  Germany  (1592-1632),  based  at  Keplar'a  Inspira¬ 
tion,  load  to  new  developaenta  of  eoaputatlonal 
aathoaatlos.  The  Alsobra  by  Burgi  was  edited  by 
Kaplar,  especially  because  It  contributed  to  tha 
eoaputatlonal  techniques.  Many  other  laportant 
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developaenta  happened  In  Prague,  especially  In 
connection  with  the  University;  nevertheless,  I 
will  not  go  into  details  except  to  emphasize  that 
thla  progress  waa  vary  closely  related  to  the 
developaent  of  applied  aathoaatlos. 

2.  Tha  period  1918-19*5 

After  World  War  X,  Czechoslovakia  was  estab¬ 
lished  as  a  deaocratlo  republic.  The  developaent 
of  ooaputatlonal  aathaaatlea  was  closely  related 
to  applications  especially  in  engineering.  Let  aa 
aantlon  as  an  axaapla  tha  fields  with  which  I  aa 
faalliar,  the  structural  Mechanics,  elasticity, 
strength  of  notarial.  One  outstanding  scientist 
in  this  direction  was  Z.  Bezant,  professor  or  tha 
Technical  university  In  Prague.  Traditionally, 
ooaputatlonal  aathods  for  tha  analysis  of  fraaa 
constructions  ware  or  great  Interest.  Essential¬ 
ly,  these  techniques  wara  related  to  tha  direct 
and  Itaratlve  aathods  for  solving  systaas  of  lin¬ 
ear  algebraic  equations.  These  usually  sophis¬ 
ticated  aathods  ware  based  rather  on  physical  and 
engineering  intuition  than  on  aatheaatlcal  theo¬ 
ries,  beosuse  at  this  tlae  aaxlaal  simplification 
was  needed  for  any  computation.  Some  of  these 
aathods  could  be  described  today  as  the  splitting 
method,  block  Iterations,  seas  as  the  method  of 
dimensional  reduction,  ate. 

Approx lasts  aathods  for  analysis  of  plates  and 
shells  based,  for  example,  on  Fourier  method,  se¬ 
ries  method,  etc.,  were  typical  for  solving  par¬ 
tial  differential  equations.  Various  solution  me¬ 
thods  had  tha  character  of  finite  differences  de¬ 
rived  on  physical  grounds  by  "spring  analyses." 
Let  us  mention  that  Cauchy's  spring  nodal  of  an 
elastic  medium  can  be  Interpreted  as  finite  dif¬ 
ference  scheme  for  Laae-Navler  equations  with 
Poisson  ratio  v  -  1/3.  Various  methods  for  solv¬ 
ing  nonlinear  problems,  eigenvalue  problems,  etc., 
were  developed  In  connection  with  buckling  and 
stability  considerations  In  general.  In  mechan¬ 
ical  engineering  various  aathods  were  developed  In 
connection  with  vibration  problems,  etc. 

The  first  mathematical  book  [1,  193*]  written 
In  193*  by  two  professors  of  mathematics  at  the 
Technical  University  in  Prague  became  a  widely 
used  text.  This  book  covered  essentials  of  numer¬ 
ical  analysis  In  a  relatively  accurate  and  de¬ 
tailed  manner.  Although  this  book  did  not  brake 
new  grounds  or  Introduced  new  approaches.  It  be¬ 
came  a  aajor  source  of  education  In  computational 
mathematics  and  in  computational  research  in  en- 


ginaering  application*  in  Czechoslovakia. 

Czechoslovakia  was  a  highly  developed  Indus¬ 
trial  oountry .  The  Skoda  enterprises,  an  indus¬ 
trial  concern,  supported  a  theoretical  department 
which  was  heavily  Involved  in  coaputations.  Thanks 
to  that,  Czechoslovakia  had  a  broad  and  firn  tra¬ 
dition  in  applied  aathenatics  and  through  it  in 
computational  set  hods. 

It  is  interesting  to  compare  the  scientific 
situation  in  Czechoalovakia  and  Poland.  Without 
any  doubts,  Poland  was  a  superpower  in  pure  mathe¬ 
matics  during  this  time;  it  was  in  the  absolute 
forefront  of  the  world  research  in  developing  such 
mathematical  fields  as  functional  Analysis,  Real 
Analysis,  Topology,  etc.  On  the  other  hand,  in  my 
opinion,  the  level  of  applied  mathematics  was 
higher  in  Czechoslovakia  than  in  Poland. 

In  the  Fall  of  1938,  Czechoalovakia  was  crip¬ 
pled  by  the  Munich  treaty;  on  March  15,  1939, 
Hitler  occupied  Bohemia  and  Moravia,  the  indus¬ 
trial  western  part  of  Czechoslovakia,  and  created 
a  puppet  state  of  Slovakia  from  the  eastern  part 
of  Czechoslovakia.  In  other  words,  Czechoslovakia 
ceased  to  exist. 

On  November  17,  1939,  Hitler  dosed  all  unl- 
veralties  to  prevent  the  higher  education  of  the 
Czech  population.  Universities  were  closed  until 
the  end  of  the  war  and  the  collapse  of  Hitler's 
Germany.  This,  of  course,  had  a  profound  effect 
in  the  development  of  science  in  general,  and 
mathematics  in  particular.  Although  there  were 
underground  seminars  and  some  mathematical  work 
and  some  more  elementary  publications  were  somehow 
published,  an  entire  generation  of  scientists  (6-8 
years  period)  was  lost.  (Same  effects  of  this 
will  be  discussed  in  the  following  sections.) 

3.  The  early  post  war  period.  Period  of  basic 

education 

Almost  immediately  after  the  end  of  the  war, 
the  Universities  were  opened  and  maximal  efforts 
started  to  fill  the  gap  created  by  the  closing  of 
the  schools  for  six  years.  Shortened  studies  were 
designed  to  fill  the  gap  as  quickly  as  possible. 
Basic  lectures  were  given  in  theaters  far  1500- 
2000  students.  This  emergency  education  had  sur¬ 
prisingly  good  effects  because  of  the  high  motiva¬ 
tion  of  the  students  and  teachers.  In  three  to 
four  years  the  major  part  of  the  educational  gap 
was  closed,  especially  In  the  education  of  engi¬ 
neers,  teachers,  medical  personnel,  etc.,  but 
could  not  and  was  not  completed  In  the  field  of 
science  and  in  the  education  of  scientists. 

In  February  19*8,  the  Communist  party  took 
over  the  government.  The  pattern  of  Soviet  organ¬ 
ization  was  applied  In  Czechoslovakia  including 
scientific  education  and  research.  Already  In 
19*9  the  institution  of  "Aspirants"  was  estab¬ 
lished.  "Asplrantura”  was  an  organization  for 
graduate  studies  In  and  outside  the  universities. 
Aspirants  were  awarded  fellowships.  Almost  at  the 
same  time,  preparations  for  the  foundation  of  the 
Academy  of  Sciences  (Soviet  style)  was  started. 

In  mathematics  the  major  responsibility  for 
the  education  of  aspirants  was  given  to  E.  tech, 
professor  at  Charles  University,  a  well  known 
topologist.  He  gathered  about  a  dozen  of  the  best 
and  moat  promising  young  students,  graduates  from 
the  universities,  and  led  their  scientific  educa¬ 
tion.  Let  me  mention  a  few  names  from  this  group 


which  became  well  known  In  mathematics  In  and  out¬ 
side  of  Czechoslovakia.  I.  Bebuika  (Numerical  and 
Applied  Mathematics),  M.  Fiedler  (Theory  of  Matri¬ 
ces),  J.  Kurzwell  (Theory  of  Ordinary  Differential 
Equations),  V.  Ptak  (Functional  Analysis),  0. 
Vejvoda  (Differential  Equations),  M.  Zlamal  (Fi¬ 
nite  Element  Method).  Under  the  leadership  of  E. 
2ech,  the  best  Czechoslovak  mathematicians  partic¬ 
ipated  in^thls  program.  I  would  like  to  mention 
especially  V.  Knichal,  V.  Kofi nek,  professors  at 
Charles  UrU  varsity  in  Prague,  F.  Vyjlchlo,  Profes¬ 
sor  at  Technical  University,  0.  Boruvka,  Professor 
at  the  University  In  Brno.  This  group  of  students 
and  their  teachers  were  a  congenial,  dedicated 
group  of  the  highest  quality.  I  have  not  seen 
afterwards  anywhere  In  the  world  such  a  congenial 
group  of  students  and  teachers. 

Professor  E.  Cech ,  although  a  pure  mathemati¬ 
cian  with  basic  Interest  In  topology  and  geometry, 
had  very  broad  views  which  he  imposed  on  the  group 
together  with  his  dedication,  hard  work  and  Inter¬ 
est  in  every  aspirant  (student).  E.  Cech  Insisted 
that  all  of  his  "aspirants"  became  familiar  with 
numerical  methods.  To  this  end,  he  obtained  from 
the  Soviet  Union  seme  old  copies  of  the  book  of 
Kantorovich  Krylov  [2,  1936],  which  was  well  known 
in  the  Soviet  Union  and  was  translated  later  In 
the  Mast.  Because  the  copying  machine  did  not 
exist  at  that  time  in  Czechoslovakia  with  the  ex¬ 
ception  of  the  ditto  sheet  machine,  E.  Cech  trans¬ 
lated  and  dictated  it  to  his  secretary,  so  that 
the  entire  book  was  typed  and  by  ditto  technology 
given  to  his  aspirants.  This  and  similar  Cech ' s 
acts  were  typical  of  his  dedication.  Neverthe¬ 
less,  it  is  necessary  to  say  that  Prof.  E.  Cech 
was  a  highly  demanding  person,  completely  "ob¬ 
sessed"  by  mathematics  (in  the  best  sense  of  the 
word)  who  permanently  challenged  his  students 
Individually  and  as  a  group  almost  In  a  dictato¬ 
rial  fashion.  In  retrospect,  one  has  to  admire 
more  and  more  his  mathematics,  dedication,  wisdom 
and  what  he  gave  to  "his”  youngsters  (with  or 
without  their  consent). 

E.  Ceoh  also  Insisted  that  the  aspirants  will 
get  basic  education  In  computer  technology  and  its 
use.  He  arranged  for  lectures  by  Prof.  A.  Svoboda. 
A.  Svoboda  worked  In  the  field  of  electronics  in 
the  United  States  during  World  War  II.  He  returned 
to  Czechoslovakia  In  19*6  and  went  back  to  the 
United  States  In  1966.  A.  Svoboda  was  the  leader 
in  the  development  of  computers  In  Czechoslovakia. 
Under  his  leadership,  a  design  and  Implementation 
of  a  unique  relay  computer  was  made  (tubes  were 
not  available  at  this  time).  Svoboda 'a  machine 
called  SAP0  was  a  triplet  machine  with  three  a- 
rlthmetlc  units  which  after  every  operation  (made 
simultaneously)  "voted"  and  the  majority  vote  was 
used  as  the  answer.  The  programming  was  a  5  ad¬ 
dress  system.  The  computer  SAP0  had  many  unique 
features.  Unfortunately  It  was  completed  when  the 
next  generation  (tubes)  was  already  In  full  swing. 

During  this  period,  work  seminars  were  rou¬ 
tine.  Teachers,  as  well  as  students,  were  In¬ 
volved  in  these  seminars.  I  remember,  for  exam¬ 
ple,  the  work  In  a  paper  by  Goldstlne,  Neisaann  [3, 
19*7]  which  convinced  us  that  there  was  no  hope 
that  elimination  method  could  and  would  be  used  In 
the  future  for  matrices  larger  than  100  (what  a 
wrong  conclusion!)  Another  paper  having  big  Im¬ 
pact  was  the  one  by  Courant,  Friedrich  and  Levy 
[*,  1927]  which  was  analyzed  In  every  detail;  E. 
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«.  Ihlldtag  the  Nathanetlcal  Institute  of  the 

Csohoalovak  loadwy  of  Solenoe 

In  the  early  firtlea,  the  Mathaaatloal  Insti¬ 
tute  of  the  Caeehosloeak  Aoadaay  was  established. 
K.  Ceoh,  V.  Knlehal,  J.  Novak,  F.  vySlchlo,  to¬ 
gether  with  aowe  of  the  prevloue  aspirants,  played 
a  prowl  neat  role  la  loading  the  Institute.  New 
res ear  ah  groups  were  built  end  another  generation 
of  youag  reseer nhare  educated. 

In  the  field  of  Applied  and  Nusarleal  He  the¬ 
se  ties  and  Partial  Differential  Equations,  X. 
Ba  busks  and  K.  hektorys*  be  pane  eery  aetlve  in 
oollaboratloa  with  Prof.  P.  VySiohlo. 

The  aaln  oaphasis  la  this  direction  of  applied 
•a  thane  ties  was  asohanlos  of  solids  and  partial 
differential  equations,  espsolally  of  slllptlo 
type.  The  aala  direotloa  was  the  relation  between 
nodarn  asset  aathenstlos  and  appliostlons  with  on - 
phssfts  on  eons  true  tlve  approaches  whloh  oould  be 
used  far  oonersto  solution  of  prahlens.  One  of 
the  result  of  this  effort  was  the  book  [i,  19533* 
The  basis  of  this  book  was  the  theory  of  analytic 
flotations  of  soapier  variables  in  the  spirit  of 
the  Musohsllahvlli  theory.  This  philosophy  of  the 
honest  no  thane  ties  la  application  later  led  to 
the  book  [5,  19663  by  K.  hektorys  and  ooworkera. 

The  above  philosophy  in  its  purest  fora,  and 
influenced  by  Nourbakl,  lad  to  sane  effort  (o.g. 
by  V.  Knlehal  and  others)  to  areata  an  axioms  tlo- 
preeise  ays ten  of  applied  aathenstlos.  This  ef¬ 
fort  did  not  aoeoapllah  too  ernoh. 

The  early  post-war  period  (I  oall  It  'period 
of  education’)  ended  roughly  In  195*  when  the 
Hathsnstloal  Institute  was  flrnly  established. 

5.  The  Project  Orllk 

The  projoot  Orllk  wea  an  isportant  nllestone 
in  the  developnent  of  oonputatlonal  and  applied 
natheaatles  In  Cseohoelovakla.  This  project  was 
nontlonod  aa  the  one  of  the  principal  aohleveaenta 
of  the  Csechoslovak  lee deny  of  Sciences  on  the  oc- 
oaslon  of  the  oolebratlon  of  30  years  of  Its  foun¬ 
dation,  and  in  the  publication  [8,  1986]  on  the 
oooaslon  of  forty  years  of  poet-war  nathcaatles. 

The  project  Orllk  wea  a  large  scale  ooaputa- 
tlonal  project  (although  atlll  perfomed  on  desk 
os leu  la  tors)  whloh  oould  be  character 1 ted  as  the 
transition  fron  the  preoenputor  to  the  oonputer 
era  In  Cseohoelovakla,  aee,  o.g.  C9,  19863.  This 
projoot  had  a  profound  la pact  and  was  character¬ 
ised  by  the  principles  whloh  after  thirty  years  are 
atlll  the  center  of  Interest  In  ecnputational  and 
applied  nathenatlcs  In  the  United  States  and  else- 
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The  research  project  Orllk  was  related  to  the 
proposed  building  of  the  largest  den  la  Csschoslo- 
vakia  located  about  *0  a  south  fron  Prague  on  the 
river  Vltava.  Tbs  dan  wea  of  oonersto  gravita¬ 
tional  type,  about  *00  ft  higi.  The  projoot  Orllk 
was  an  lntagrated  couple*  reseerch  in  wa thane ties, 
engineering  and  notarial  aolonoe  Consent,  con¬ 
crete)  .  The  leader  of  the  nathonatloal  part  was 
I.  NebuXka.  of  the  engineering  part  Prof.  L. 
Hejslfk  (Professor  of  Tech.  University  Brno),  and 
of  the  toohnologtoal  part.  Dr.  J.  Jlrsak.  The 
projoot  was  a  teas  work  and  Included  a  large  staff 
of  people  working  on  desk  calculators. 

The  aaln  technical  prot) Ion  was  that  the  con¬ 
crete  releases  a  significant  anount  of  heat  dur¬ 
ing  hardening.  Sinul tenuously,  the  hardening, 
whloh  deponds  strongly  on  the  toaperature,  changes 
significantly  the  notarial  properties,  e.g.  elas¬ 
ticity  nodule,  ereap  and  relaxation  properties, 
etc.  This  leads  to  the  creation  of  significant 
stress  state  whloh  la  "frozen  In"  during  the  hard¬ 
ening  and  later  oould  lead  to  dangerous  and  seri¬ 
ous  eraoks.  The  effects  of  this  type  oould  be 
oontrollad  by  a  proper  technology  of  building  and 
of  notorial  proportion.  The  large  dans  In  the 
United  States  used  a  cooling  ays ten  by  pipes  In¬ 
serted  In  the  dans.  The'  baslo  questions  of  the 
rssesroh  worst  a)  Wist  are  the  effeeta  of  various 
building  proosdvros  on  the  possible  cracks?  Is  It 
necessary  to  use  pipe  cooling,  etc.?  Could  the 
-t  orooks,  if  any,  be  expected?  b)  How  the  proper¬ 
ties  of  the  oonorote  lnfluenoe  the  undesired  ef¬ 
fects  of  building.  1st or  functions  of  the  dan, 
etc.?  Based  on  the  rosoaroh  results,  the  dan  was 
built  without  cooling  by  a  relatively  quick  build¬ 
ing  sohedule  In  blocks  about  12  ft  high.  The  dan 
behaved  as  predicted  and  serves  well  Its  purpose. 
Neoulta  of  the  analysis  were  presented  at  the  dan 
world  congress  in  1958,  and  were  Included  dlstln- 
gulshly  In  the  congress  reporter's  report.  Sane 
teehnlasl  conclusions  are,  e.g.  In  [10,  1958], 
[11,  19613. 

The  essential  novelty  was  the  oaphasis  on 
Integrated  approach  and  the  reliability  of  the 
conclusions.  The  reliability  aspects  were  divided 
In  the  following  groupss 

a.  reliability  of  nathonatloal  nodel, 

b.  reliability  of  available  Input  data, 

o.  reliability  of  the  nunerioal  aethod  and 
principles  of  its  selections, 

d.  reliability  of  the  arlthaetlc  oonputatlons 
(round  offs)  (bean use  nlninal  nusber  of 
used  digits  were  essential  for  computa¬ 
tions  on  desk  calculators). 

These  questions  were  directly  and  Indirectly  ad¬ 
dressed  in  a  series  of  theoretical  and  engineering 
papers  and  reports. 

The  proto  Isa  was  highly  nonlinear  and  three- 
dlaenslonal.  Beoause  three-dinenalonal  solution 
was  out  of  the  cpiestlon  for  obvious  reasons,  a 
series  of  two-dimensional  problens  were  solved  and 
combined  approxiaately  Into  three-dimensional  ones 
by  a  sort  of  splitting  up  approach.  Let  me  ex¬ 
plain  now  some  or  the  problens  (In  a  simplified 
way). 


K. 


Nektorya  Is  the  author  of  [53  and  [73,  which  are  well  known  In  the  United  States. 
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Bar*  u  is  th*  temperature,  F  th*  Intensity  of 
the  host  or*at*d  by  hydration,  4  a  fiotlv*  tin* 
(as*)  ia  which  th*  saa*  aaount  of  hast  was  pro¬ 
duced  aa  wh*n  th*  t«np«ratur*  would  b*  fixed 
(about  70*P).  This  fiotlv*  tin*  characterizes  the 
state  of  the  ohanioal  reaction.  Tha  coefficient* 
e(u,«)  and  a(u,p)  war*  found  so  alldly  depend¬ 
ant  on  u,«  that  average  values  war*  us*d.  Th* 
characterisation  of  F(u,4)  was  essential.  A 
special  oar*  was  devoted  to  th*  laboratory  experl- 
nanta.  Finally,  the  above  aantloned  aodel,  based 
an  a  ohanioal  nodal  of  hydration,  was  aooepted  and 
a  differential  aquation  (1a,b)  was  designed  and 
used.  Th*  date  w ere  obtained  by  the  a*aauran*nt  of 
tha  heat  release  in  the  period  (0,t)  under  con¬ 
stant  tenperatures  and  in  an  adiabatic  state.  The 
ooaputetlan  of  the  1 nor e* ante  In  F  was  organised 
ao  that  th*  total  hast  was  exaotly  preserved.  This 
was  very  essential  for  the  reliability. 

The  technology  of  the  building  consisted  in 
quick  production  of  blocks  about  12  ft  high  with 
tin*  intervals  T  in  between.  The  ache**  ia 
shown  in  Figure  1.  To  slapliry  the  problen,  a 
periodic  solution  (in  tins  and  apace)  waa  ana¬ 
lysed.  It  has  been  shown  that  the  solution 
quickly  approaches  th*  state  u(t+T,x,y«d)  - 
u(t,x,y),  Oi  t  1  r  and  this  state  was  nuneri- 
oally  oonputed  [12,  I960]. 

The  nunerloal  netted  waa  essentially  th*  fi¬ 
nite  difference  netted  with  the  schene  derived  by 
the  oell  Integration  identity  principle*  guaran¬ 
teeing  th*  balano*  condition.  This  technique  was 
dose  to  the  teohnlque  of  Marchuk's  identity,  ela¬ 
borated  later  in  [13,  1966]. 

An  essential  feature  which  was  introduced  nuoh 
later  in  the  finite  elenent  netted  under  the  nan* 


'special  el  seen  ts*  was  uaed  in  th*  eoaputatlons. 
In  th*  presence  of  ooollng  pipes  there  waa  a  sig¬ 
nificant  teat  sink.  Hence,  the  solution  waa  writ¬ 
ten  in  the  for* 


u(x,y,t)  -  v(x,y,t)  ♦  w(x,y,t) 


wher*  w(x,y,t)  waa  the  linear  solution  of  a 
point  source  (nor*  precisely  single  olrol*  source) 
with  th*  Intensity  o(t)  (which  waa  the  oonputed 
intensity  of  ooollng).  Ftaetlon  v  waa  deter- 
ained  by  finite  difference  netted  aa  explained 
above  and  the  hydration  heat  was  Included  In  this 
ten*.  (For  th*  stationary  solution  exactly  th* 
wetted  of  special  elenents  was  obtained.) 


2)  The  free zlna  problen.  Th*  building  of  th* 
dan  had  to  oontlnu*  during  winter  when  freezing  of 
the  oonorete  in  th*  beginning  phase  of  hardening 
oould  create  a  serious  da  sage.  At  aost  the  con¬ 
crete  la  allowed  to  freeze  for  a  abort  tine  at  a 
depth  of  1  -  2  in.  Th*  wooden  siding  for  laying 
the  oancret*  serves  also  as  insulation,  and  the 
freezing  ooours  when  the  siding  Is  noved  in  the 
next  building  cycle.  Th*  aaln  approach  was  here 


to  solve  a  stochastic  problen  for  Equation  (1 ) 
when  th*  boundary  oondltlons  are  a  stochastic 
function  -  the  outside  tenperatur*.  The  sain  val¬ 
ue  and  dispersion  for  the  desired  lnforaation  were 
oonputed.  The  theoretical  base  was  described  in 
[It,  1961].  (Let  us  renark  that  today  a  large  re¬ 
search  project,  sponsored  by  RASA  Lewis,  solving 
this  problen  with  stochastic  input  data  is  in  pro- 
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Fig.  1.  Schenatic  state  of  progressed  dan. 

3)  The  elasticity  problen.  Given  the  tenper¬ 
atur*,  the  thernos tresses  were  cow pu ted.  The  es¬ 
sential  problen  was  the  fornulatlon  of  the  problen 
with  respect  to  naterlal  properties  Including 
change  of  elasticity  nodule  and  creep  (relaxation) 
properties,  etc.  A  rheological  aodel  based  on  a 
description  of  the  chenlcal  process  of  hardening 
was  designed  and  tested  In  the  laboratory. 

Th*  nuaerlcsl  solution  was  based  on  a  series 
of  plane  problems  in  the  spirit  of  splitting  up 
nethods.  In  this  phase,  J.  le&s  contributed  sig¬ 
nificantly  to  this  research.  Anong  others,  the 
theoretical  papers  [15,  1958]  [16,  1959]  are  dl- 
reotly  related  to  this  work.  The  nonograph  [17, 
1967]  by  J.  MeSas  Is  the  only  basic  nonograph 
which  does  not  avoid  unsnooth  dona 1 ns .  This 
nonograph  and  other  results  or  J.  HeSas  are  well 
known  in  the  West.  Various  Iterative  nethods  were 
used  in  the  connection  of  splitting  the  problen 
Into  two  dlnenslonal  ones.  Let  us  sent  Ion  one  of 
the,  type  of  Schwarz  alternating  algorlthn. 
Mathenatloally,  the  aaln  generalization  used  was 
based  on  the  following  functional  analytical  frane 
(which  Is  today  nore  or  less  standard),  formulated 
her*  in  the  slnplest  forw: 

Let  ,  ?2  be  projection  operators  on  the 
subspeces  .Sj  c  h.  Then  (P2?2)n  converges 
polntwise  to  the  projection  onto  S1  n  S^. 

A)  Error  control .  The  basic  Idea  of  the 
error  control  of  the  numerical  netted  was  to  In¬ 
terpret  the  nunerical  solution  as.  exact  solution 
of  a  problen  with  slightly  different  Input  data. 
The  natheaatieal  models  were  verified  by  computa¬ 
tion  of  son*  slnple  laboratory  experiments.  The 
round-off  error  was  analyzed  In  a  way  close  to 
that  explained  later  In  the  monograph  [13,  1968] 
by  a-sequences.  In  the  project  Orllk,  a  team  of 
researchers  was  Involved.  In  addition  to  those 
already  mentioned,  I.  BabuSka,  L.  Mejzllk, 
J.  Jirsak,  E.  vitasek,  J.  NeSas,  other  researchers 
participated,  especially  K.  Rektorys,  M.  Prager, 


f.  vySloblo.  Various  publications  and  reports, 
v  union  dirootly  or  lndireetly  wars  rolatod  to  the 
projaet,  wore  published  during  this  tine. 

6.  Tha  raaaaron  in  the  optimisation  of  the  numer- 

ioal  nathods,  nunarloal  stability  and  numeri¬ 
cal  methods  in  general 

Ouring  the  sixties  (1964,  1967),  conferences 
devoted  to  nunarloal  mathematics  mere  organised. 
Emphasis  urns  placed  on  the  questions  of  optimality 
of  the  selection  of  nunarloal  nethod  and  nimierical 
stability.  These  conferences,  which  took  place  in 
the  castle  Llbllce,  were  held  in  a  very  informal 
working  atmosphere.  Leading  numerical  analysts 
and  mathematicians  from  east  and  west  partici¬ 
pated.  Let  me  mention,  among  others,  M.  S. 
Bachvalov,  C.  Golub,  P.  Henri ci ,  G.  I.  Marchuk, 
F.  Olver,  S.  L.  Sobolev,  A.  N.  Tlchonov.  These 
conferences  were,  in  my  opinion,  the  very  first 
meetings  in  the  world  concentrating  specifically 
on  the  questions  of  optimal  selection  of  the  nu¬ 
merical  method.  The  various  aspects  of  optimal¬ 
ity,  theoretical  and  computational  were  discussed. 
Some  Ideas  and  results  related  to  this  direction 
obtained  In  Czechoslovakia  were,  for  example,  pre¬ 
sented  In  Cl 3*  1966]. 

B.  Computational  mathematics  1m  USSR 

In  this  section  I  will  make*  a  few  subjective 
comments  about  the  development  of  computational 
mathematics  In  USSR  up  to  the  mid  1990.  For  a 
systematic  survey,  we  refer  to  [18.  19*8]  and  [19, 
1959]. 

The  theory  of  approximate  methods  has  a  long 
tradition.  For  example,  the  idea  of  the  Galerkin 
method  was  introduced  in  1915  In  [20,  1915].  The 
Rltz  method  was  investigated  in  a  series  of  papers 
of  Krylov  and  BogolJUbov.  See  e.g.,  [21,  1917] 
[22.  1917],  (23,  1927],  (2*.  1931].  The  Galerkin 
method  was  investigated  by  various  authors  in  the 
pre-war  period.  The  book  of  Kantorovich  and 
Krylov  [2,  1936]  is  likely  the  first  comprehensive 
book  about  the  mawrloal  solution  of  partial  dif¬ 
ferential  equations.  (After  the  war  this  book  was 
translated  into  many  languages.) 

The  Faddejeva's  sonography  [25,  1950]  Is  like¬ 
ly  the  first  comprehensive  book  about  the  methods 
of  linear  algebra.  (It  was  later  translated  into 
English. )  Michlln's  work  and  books  (e.g.  [26, 
1950s  27,  1952]  and  others)  about  the  variational 
methods  were  important  contributions  to  the  theory 
of  variational  methods  and  computational  ap- 
proachea. 

1 .  Variational  methods 

As  I  have  already  mentioned,  the  variational 
methods  were  investigated  by  many  authors.  The 
investigations  addressed  both  the  Rltz  nethod 
baaed  on  a  minimization  of  a  quadratic  functional 
as  well  as  the  Galerkin  method  (sometimes  called 
methods  of  moments  or  weighted  residuals)  with  the 
sane  or  different  trial  and  test  spaces.  The  re¬ 


sults  related  to  the  applications  of  a  minimiza¬ 
tion  are  using  the  Friedrichs  extension  of  the  op¬ 
erator  to  a  self  ad  joint  one.  This  direction  was 
utilised  by  Michlln  in  many  of  his  papers  and 
books,  and  Michlln  was  likely  the  first  who  used 
the  tern  "energy  space.”  An  important  role  played 
the  analysis  of  the  energy  space  and  the  question 
to  what  Sobolev  space  (in  today's  terminology)  it 
is  equivalent.  For  example,  in  [27,  1952]  this 
question  is  analyzed  for  basic  problems  of  the 
elasticity  theory.  For  the  nixed  problem  (e.g. 
free  friction  contact  boundary  condition)  the 
equivalency  was  analyzed,  e.g.  in  [28,  1951].  The 
characterization  of  the  energy  space  for  Poisson 
problem  on  an  infinite  domain  was  discussed  in 
[29,  1953].  The  convergence  of  the  Rltz  nethod  in 
the  energy  space  is  then  directly  related  to  the 
best  approximation.  An  effort  was  made  to  analyze 
the  convergence  in  the  stronger  norms  |u|  - 
(Au.Au)*  (see,  e.g.  [30,  1956])  or  weaker  norms 
as  Ml.  (see,  e.g.  [31.  19*1].  The  convergence 
of  the  Treftz  method  was  analyzed  in  detail  In 
[26,  1950]. 

The  Galerkin  nethod  and  general  method  of  mo¬ 
ments  (also  with  different  trial  and  test  func¬ 
tions)  for  Integral  equations  were  studied  in  many 
papers  by  Krylov  and  his  coworkers.  See,  e.g. 
[23,  1927],  [32,  1931].  In  applications  to  dif¬ 
ferential  equation,  Petrov  [33,  19*0]  used  the 
different  trial  and  test  spaces,  and  the  term 
Galerkln-Petrov  nethod  is  used  sometimes  today. 
Keldys  [3*.  19*2]  applied  this  nethod  to  a  non- 
sal  fad  Joint  boundary  value  problem  for  ordinary 
differential  equation;  this  paper  very  likely  was 
the  first  one  establishing  the  convergence  of  the 
nethod  in  general  setting  when  applied  to  an  spe¬ 
cific  problem.  The  convergence  or  the  Galerkin 
nethod  was  established  by  Michlln  tor  the  opera¬ 
tors  of  the  form  A  -  Aq  ♦  K  where  Aq  is  posi¬ 
tive  definite  selfadjolnt  operator,  and  A q’k  is 
compact  in  the  norm  (Aqx.x)V  See  [35,  19*8], 

[36,  1950],  [37,  1957].  In  [38,  19*8],  a  general 

functional  analytic  scheme  of  numerical  nethod  was 
discussed  by  Kantorovich.  See  e.g.,  [39,  I960]. 
The  main  idea  is  roughly  the  following.  Let  us  be 
Interested  in  Kx  •  y  with  x  €  X,  y  €  Y.  Then 
the  numerical  nethod  solves  essentially  K),xh  - 
yh  where  h  Is  a  parameter,  h  ■»  0  and  xh  €  X 
y.  €  f.  There  is  a  one  to  one  napping  <pj,  of 

X  onto  X  €  X  and  p),  of  ?  onto  fey. 

Then- one- would  like  to  achieve  that  is 

close  to  the  solution  of  the  original  problem.  For 
that,  one  has  to  essentially  achieve  that  phK  - 
KtfPfl  is  small.  In  [38,  19*8]  this  approach  was 
applied  to  a  large  class  of  illustrative  problems. 

Collocation  nethod  obviously  can  also  be  un¬ 
derstood  as  nethod  of  moments  and  has  been 
treated,  e.g.  in  [39,  1960],  In  the  frame  of  the 
above  mentioned  approach.  A  nethod  which  Is  very 
close  to  the  collocation  was  applied  in  [*0, 
195*],  [*1 ,  1956]  by  Vlshlk.  In  an  abstract  form, 
the  Galerkin  nethod  and  nonlinear  problems  and  a 
discussion  of  the  approximate  nethod  are  given  by 
KrasnoselsklJ  in  [*2,  195*]  and  In  some  of  his 
other  papers. 


*1  give  here  the  references  to  the  originals  In  Russian.  Translations  of  many  of  these  papers  and  books 
are  now  available. 
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*.  Halt*  aAM 

Tin  tMi«  ttwry  of  tin  finite  differ  no#  no- 
Uni  eepeolmlly  rotated  to  the  stability  is  in  Uw 
book  by  ajnbonkjj  Filippov  [ij,  1996].  A  handbook 
of  finite  difference  aeh— aa  far  partial  differen¬ 
tial  equations  was  nr  it  ton  by  Panov  CM.  1951]. 
For  hyper belie  equatlona  there  la  a  aeries  of  re¬ 
sults  of  tedysoMdcaJa  and  her  ooworfcera.  See,  e.g. 
CbS.  1952],  £t6,  199*].  CM.  1953]. 

In  the  ease  of  elliptic  equations,  early  works 
are  given,  for  avenple,  in  CM,  195*].  C*9,  19*7]. 
Per  applications  of  finite  differ onoe  for  parabol¬ 
ic  equation,  we  refer,  for  evanple.  to  the  work  by 
Kamynin  CSO,  1953]. 

The  general  eigenvalue  tree  font  by  the  finite 
difference  uethod  is  given,  for  exaaple,  in  C51. 
195*]. 

3.  Nunerical  treatment  of  differential  equatlona 

In  the  previous  sections  sons  early  works  were 
presented.  They  played  (by  the  subjective  Judge- 
sent  of  the  author)  iaportant  roles  in  the  devel- 
opaent  of  the  theory  of  the  iwerlcal  uethod.  It 
is  interesting  to  aentlon  that  although  the  theory 
of  variational  uethod  was  eery  advanced ,  the  entire 
direction  of  the  finite  eleuent  uethod  was  for  a 
long  tine  neglected,  and  the  aaln  euphasla  was 
placed  on  finite  difference  uethod.  It  seeas  to 
be  characteristic  the  finite  eleuent  uethod  was 
aalled  until  recently  variational  rinlta  dif¬ 
ference  uethod. 

Finite  difference  uethod  was  later  analysed  in 
the  works  of  SanarsklJ,  Oodunov  and  uany  others, 
and  uany  uonographa  and  test  books  are  available 
today.  In  these  works  the  euphasla  is  placed  on 
the  theory.  The  discussions  of  computational  as¬ 
pects.  nioMrleal  exper  loan  tat  Ion.  analyaes  of  the 
performance  of  the  uethod  on  banchaark  problems 
are  very  rare.  Very  likely  this  situation  is  re¬ 
lated  to  the  state  of  the  oouputer  technology  In 
USSR,  nevertheless,  the  ccuputer  situation  stim¬ 
ulated  varloua  special  methodologies  as  splitting 
up  uethoda.  and  various  "tricky"  Iterative  proce¬ 
dures  which  were  used  in  scientific  computations. 
In  the  area  of  uatheuatlcal  modeling  and  scien¬ 
tific  computations,  iaportant  works  have  been  done 
by  0.  I.  Marchuk  and  his  coworkers  in  uany  papers. 
The  first  of  his  books  (see  [52,  1958]  is  address¬ 
ing  modeling  and  computational  uethoda  in  reader’ 
analysis.  It  is  interesting  to  aentlon  that  the 
Idea  of  preconditioning— credited  to  Buljaev — is 
nontloned  there. 

I  only  sectioned  very  few  papers  and  results; 
nevertheless,  hopefully,  they  give  soue  illustra¬ 
tive  picture  of  the  character  of  the  research  In 
the  USSR  in  the  early  post-war  period. 

Cl]  tasks,  V.  Hruska ,  v..  Theory  and  practice  of 
oouputations,  Prague,  193*  (in  Czech). 

[2]  Kantorovich,  L.  V.,  Krylov,  V.  I.,  Approxi¬ 
mate  solutions  of  partial  differential  equa¬ 
tions,  Moscow,  Leningrad,  1936,  pp.  1-588 
(In  Russian). 

[3]  von  Neumann ,  4,  Goldstlne,  H. ,  Nunerical 

Inverting  of  matrices  of  high  order,  Bull 
tear.  Math.  Soo.  53  (19*7),  pp.  1021-1099. 

(*]  Courant,  R.  Friedrichs,  K.  0.,  Levy,  H. , 


Uber  die  pertiellen  Differenzenglelchungern 
der  Physik,  Hath.  Ann.  100,  1928-1929.  pp. 
3*-7*. 

[5]  Rektorys,  K.,  Survey  of  Applicable  Mathemat¬ 
ics,  Prague,  1966. 

C6]  Bmbuika,  I.,  Rektorys.  K.  Vyelchlo,  F., 
Ma theme tlsche  Elastizltitstheorie  der  ebenen 
Problene,  Academia- Ver lag.  Barling,  I960. 

C7]  Rektorys,  K.,  Variational  methods  In  mathe- 
uatlcs,  science  and  engineering,  Dordrlcht, 
Boston,  D.  Reldel  Pub.  Co.,  1977. 

[8]  Development  of  the  mathematics  In  Czecho¬ 
slovakia  in  the  period  1 9*5-1 985  and  further 
perspectives,  Charles  University,  Prague, 

1966.  1-217  (In  Czech). 

£9]  Marek,  I. ,  Approximate  and  numerical  methods 
In  [8]  (1985),  pp.  127-1*3  (In  Czech). 

[10]  Bebuska ,  I.,  Mejzlfk,  L. ,  Calculation  and 

measurement  of  termal  stresses  In  gravity 
dams,  VI  Congress  des  Grandes  Barrages.  New 
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